Dual proximity effect near superconductor-insulator transition 
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We show that quantum vortex-loop proliferation may be one possible explanation of the super- 
long-range proximity effect observed in an insulating underdoped cuprate, using the dual theory of 
quantum vortices. As a test of this scenario, we propose that a dual proximity effect experiment 
can confirm the superfluid motion of the quantum vortices in the vortex-proliferated insulator and 
can measure the divergent correlation lengths near the superconductor-insulator transition. 
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Cuprate superconducting materials exhibit a rich 
structure of phases depending on the charge-carrier dop- 
ing and the temperature. Most notably, they exhibit su- 
perconducting, antiferromagnetic (AF) Mott-insulating, 
and T > T c strange metallic phases. To this date, how- 
ever, the nature of the intervening phases between the 
superconducting state and the AF Mott-insulating state 
is not understood [fjj . In attempts to understand this re- 
gion of the phase diagram, there have been a number 
of theoretical proposals on the mechanism of destruc- 
tion of superconductivity such as the superconductor- 
insulator transition into the quantum-disordered insula- 
tor with vortex-loop proliferation [||-^), superconductor 
into AF insulator transition with the SO(5)-symmetric 
theory transition into novel ground states which are 
possible in two-dimensional (2D) doped Mott insulators 
or striped phases ||. Therefore, it would be desir- 
able to devise experiments which can shed light on the 
nature of the underdoped insulating state. 

In an experiment performed by Decca et al. 
Josephson junctions of superconductor- insulator- 
superconductor (SIS) type were made by generating 
photo-induced superconducting electrodes on insulating 
underdoped thin films of Yf^CusOg+a; family. Surpris- 
ingly, they observed the Josephson effect with an invari- 
ant product of the Josephson critical current and the nor- 
mal states junction resistance for separations as large as 
100 nm between superconducting electrodes, much larger 
than the superconducting coherence length £o ~ 1 nm - 
This finding indicates enhancement of the effective corre- 
lation length (£ c ) of the superconducting order parameter 
in the insulator, which can be explained by the proxim- 
ity to a continuous quantum phase transition. For in- 
stance, the SO (5) theory is one possible scenario which 
predicts such an effect [To| ]. In this paper, as an al- 
ternative scenario, we will demonstrate that the vortex 
proliferation does provide a long correlation length scale 
near SIT, which may explain the super-long-range prox- 
imity effect in the underdoped insulator. In fact, exper- 
iments on the optical conductivity ]Tl| and the thermal 
Nernst effects jl2| provided evidences for strong ther- 
mal fluctuations of unbound vortices over a wide range 




FIG. 1. A dual Josephson junction experiment of insula- 
tor-superconductor-insulator type. 



of temperatures above T c in the underdoped pseudogap 
phase, which suggest a Kosterlitz-Thouless like transi- 
tion with a broad phase fluctuation regime [ fl3|] . In ad- 
dition, an experiment on current-voltage characteristic 
of superconducting Bi2212 shows indications of a large 
density of quantum vortex pairs well below T c and far 
away from SIT E3]. Therefore, given the indications of 
strong quantum/thermal phase fluctuations in cuprates, 
vortex-proliferation is a prospective candidate scenario of 
SIT. 

We also propose an experiment which can detect a 
qualitatively distinct effect of Josephson tunneling of su- 
perfluid vortices in an insulator-superconductor-insulator 
(ISI) junction. If the SIT is induced by the vortex 
proliferation, we expect that the divergence in correla- 
tion lengths is symmetric between the superconducting 
and insulating side of the transition due to the electric- 
magnetic duality jl5| . More precisely, the vortex prolif- 
erated insulator can be considered as the bose-condensed 
state of vortices, and consequently superfluid motion of 
vortices is possible. Therefore, by choosing a supercon- 
ducting junction sufficiently close to SIT, it will be pos- 
sible to detect the super-long-range proximity effect of 
the vortex-condensate superfluid in the dual Josephson 
(ISI) junction experiment Jl6|] . In Fig. [l] we show an 
ISI junction which can be prepared by the same exper- 
imental technique used in Ref. ||. We consider apply- 
ing an external magnetic field gradient across the junc- 
tion and subsequently letting it relax, which will impose 
a time-dependent phase difference (Atp) in the vortex- 
condensate order parameter ($) across the junction. The 
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phase difference induces a dual Josephson current of mag- 
netic vortices {J) across the superconducting junction. 
The vortex current leads to a transverse electric volt- 
age difference (AVt r ) across the superconducting junction 
which can be subsequently measured. As in the super- 
conducting proximity effect experiment, we can imagine 
varying the length (5) of the junction to probe the corre- 
lation length of the magnetic vortex field $s m the super- 
conducting junction. If the superconductor is in proxim- 
ity to the SIT, we expect a large vortex correlation length 
(£ v ), analogous to the super-long-range proximity effect 
in the insulator observed in Ref. || . We show that a sim- 
ple oscillatory solution to the equation of dual Josephson 
currents can be obtained under certain conditions. 

In this paper, we use the dual theory of vortices for 
a Ginzburg-Landau description of vortex condensates in 
the insulating state. We first review the dual formulation 
P,js| ]r?| , |l8[| . Throughout the paper we use the unit h = 1 , 
and denote the space-time three- vectors as x — (x, t). 
We begin with the effective phase-only action of a two- 
dimensional (2D) zero-temperature superconductor near 
SIT coupled to the electromagnetic field A^: 

S[4>, A ^ = \f dx dt {p> t( 2 ^° + d ^) 2 l c l ~ ( 2A + W) 2 ] 
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where n s is the superconducting fcrmion density, p s is 
the superfluid phase stiffness defined as p s = n s /Amf , mf 
is the effective fermion mass, c s is the superfluid phonon 
velocity, and d is the thickness of the superconducting 
film. In a BCS superconductor, c s is related to the Fermi 
velocity vf by c s — vp/^/2. Here we chose the Coulomb 
gauge V • A = 0. 

Since we are interested in the effective theory of mag- 
netic vortices, we first separate d^ip into d^(j) = d^O + A^ 
where 9 is the spin-wave type Goldstone fluctuation and 
Afj, gives the topologically non-trivial phase gradients 
generated by magnetic vortices, and integrate out both 
the 9 and A^ fields to obtain the final effective theory of 
the A^. The A^ fields are related to the magnetic vor- 
tex current density (J?) by J* = e^ x d v A\. Then we 
transform the dynamics of single-particle paths of mag- 
netic vortices into the field dynamics by introducing a 
quantum vortex field $ which is coupled to the super- 
fluid phonon field represented by a gauge field G p . The 
resulting well-known dual form ]5|Jl7| of the effective the- 
ory of vortices follows: 



S[*,$*,G] = / d 2 x dt [\(dt + 2niG )&\* 
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where S[G^] is the action of the new U(l) gauge field 
G M . Here tC = p s k 2 /(k 2 c 2 - uj 2 + 8ire 2 p s \k\) and K T ~ 
/9 s k 2 /(A~ 2 + k 2 ) where A is the penetration depth. We 
also included the repulsive potential (it > 0) between 
vortices. The magnetic vortex rest mass is obtained from 
the energy of a pancake vortex and given by m v v 2 »s 
Ps'k In k where k = AA C , where we assumed a vortex core 
size of £o which provides the momentum (short-range) 
cutoff at A c = 7r/£o and the frequency (short-time) cutoff 
at wfA c . Here we take v v ss vf on physical grounds p8|. 
From Eq. (^), the vortex number density (n v ) and the 
vortex current can be expressed by <& as 

n v = Jo/27r = Im2$*a $ + 47rG |$| 2 (4) 
J/2n = u 2 (Im 2$*V$ + 4ttG|$| 2 ) . (5) 

In addition to Eq. (||), there is a term that couples 
the vortex currents to the background superfluid den- 
sity which acts as a dual magnetic field for vortices, but 
we assume that this does not affect the superfluidity of 
the quantum vortices in the insulating state where the 
Cooper pairs are localized. 

We now discuss the superconducting proximity effects 
in the insulating state. The action in Eq. (||), if t is ana- 
lytically continued to the imaginary time, can be viewed 
as the effective 3D Ginzburg-Landau functional of $. 
Then SIT occurs upon Bose-condensation of the <3? fields, 
when m 2 < so that ($) ^ and the correlation func- 
tion in the mean-field approximation is ($(x)$* {y)) ~ 
|$o| 2 > 0, independent of \x — y\. The correlation func- 
tion ($(x)$*(y)) can be considered as the expectation 
value of the number of magnetic vortex paths that con- 
nect x and y, and therefore ($(x)$*(y)) « |<I>o| 2 > 
implies a finite probability of arbitrarily long 2+1D vor- 
tex loops. Here we demonstrate that the Bose-condensed 
$ field introduces a new length scale related to the ob- 
served super-long correlation length . 

We first derive an effective Ginzburg-Landau func- 
tional of the superconducting order parameter (\P) in 
the vortex-proliferated insulating state. The derivation 
is straightforward once we identify the trajectories of 
Cooper pairs with dual (electric) vortex lines in the dual 
theory fl7j| . The symmetry between the Cooper pairs and 
vortices is clear from the form of Eqs. (|lj) and (|J) which 
have similar functional forms in p s e % ^ and $ respectively, 
except for the topological term dt4>n s /2 in Eq. (|]) which 
is absent in Eq. (||) and action of the gauge fields. Be- 
sides, just as a Cooper pair picks up a phase factor of 
2-7T upon encircling a vortex, a vortex picks up the same 
phase factor upon encircling a Cooper pair. Therefore, 
in the superfluid state of vortices (insulating state), a 
Cooper pair is the topological defect, playing the role of 
a vortex. Thus, we can perform a dual transformation on 
Eq. (0) to obtain an effective theory of superconducting 
order parameter, in the same way as the dual theory of 
vortices in Eqs. (0) and (0) is derived from the original 
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superconducting theory in Eq. ([|). The result is follow- 



ing: 
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Here we estimate that m c v 2 = 47r|$o| 2-y v m (-^/Co) an d 
take v c ~ v v on physical grounds, where L is the linear 
dimension of the sample. Therefore, we may take the cor- 
relation length of the superconducting order parameter 



47r|<I>o| 2 Wv ln(£/£ ), and so £ c can be 



from £~ « m c v c 

a very large number near SIT where |<I>o| 2 is small. This 
can explain the anomalous proximity effect observed in 
Ref. H , assuming that the insulating compound is in the 
vortex-proliferated state. The vortex-proliferated insu- 
lating state has a physically distinct property of vortex 
superfluidity as shown below. 

Next we consider the dual proximity effect near SIT. 
We consider an ISI junction surrounded by a supercon- 
ducting region so that magnetic flux can be trapped in- 
side. We assume that the vortex condensate density in Ii 
and I2 is l^ol 2 - We imagine applying a magnetic flux den- 
sity difference (AB) between Ij and I2 regions of the ISI 
junction and relaxing it at a certain time (t = 0) as shown 
in Fig. |], and consider the time evolution of the mag- 
netic field difference and the vortex supercurrent across 
the S region. Here the difference in the total number of 
magnetic vortices is J d 2 x(AJ ) /2n — ABL 2 /^ where 
4>a = hc/2e is a quantum of magnetic flux in a supercon- 
ductor and L is the linear size of the Ii and I2 regions. 
This condition can be implemented by adding a Lagrange 
multiplier J dt ,u v (i) / d 2 x[A J (x, t) - 27rAi?/0o] to Eq. 
(^). The chemical potential difference /i v can be ab- 
sorbed into the dual electro-chemical potential difference 
Go which is the sum of the chemical potential and the 
dual scalar potential generated by the presence of extra 
vortices introduced by the external magnetic field. Ignor- 
ing the spatial fluctuations in Go, we can absorb Go into 
the phase difference of the field $ by the transformation 
$ -> $ e iAv where 



d t A<p(t) = -2ttGo(*) 
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The difference in the vortex number density can be re- 
lated to Gq from Eq. (0) as follows: 



An v = (AJ /2ir) n 4nG (t)\$\ 



(9) 



neglecting the time derivatives of the condensate <&. The 
value of Go is determined by the condition An v = 
AB/(f>o w 4ttGo|$| 2 , which gives G ~ AB/(4w\^ Q \ 2 ^ ). 
Inside the superconducting junction (S) of length 5 and 



width L, the $s held is phase-incoherent and does not 
have non-zero vacuum expectation value. From Eq. (||), 
ignoring non-linear effects, the semi-classical $5 field 
minimizes the following action: 
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with a correlation length of £ v ~ l/m v v v . If the field 
$s ( x j t) is solved in the S region, we can obtain the vortex 
number current [I v = L^J x {x, t)/2ir] between Ii and I2 (in 
the x direction) flowing through the S region. Then the 
rate of change of the vortex number density difference 
can be related to the vortex current as follows: 



d t An v = J x (-K,t)/2TrL 



(11) 



Below we show that if we assume a nearly one- 
dimensional geometry of the junction, a simple solution 
to Eqs. (||), (||), (|I(i|), and ( |TT| ) can be obtained under 
following conditions: (a) The magnetic field difference 
(AB) is small so that |<E>| 2 in the insulators can be ap- 
proximately taken as a constant as a function of An v . 
(b) Inside the superconducting junction region, the fol- 
lowing conditions hold: I'&s/'Cvl 2 <C Id^sl 2 , which is 
valid when 6 < £ v , and \d t <&$/v F \ 2 < l^^sl 2 - In this 



case, from Eq. (10), the equation of motion for <&g(x,£) is 
simply reduced to a static Laplace equation: 9 2 $ s = 0. 
(c) The phase difference remains small (\A<p\ -C 1) due 
to the small magnitude of the magnetic field difference 
(AB). From condition (b), along with the boundary con- 
ditions $s(0) = §oe' lAv and $s(#) = $0, we can obtain 
the solution of the Laplace equation in the region S p9[ : 



§ [e^(l~x/8) + (x/5)] 



(12) 



Then one can show that there is a magnetic vortex su- 
percurrent through the region S as follows: 

J x a 2Trv 2 i(<i> s d x ^* s - $* s d x $ s ) « -£J v sinA<^ , (13) 

where J v = 47rw 2 |<i>o| 2 /(5. If 5 exceeds the correlation 
length £ v , J v will be exponentially suppressed in 5/^ v 
p9[ . From Eqs. (@) and @, and substituting the vor- 
tex current obtained from Eq. ( |l2|) into Eq. (11), we 
obtain the equation d 2 Aip = — a; 2 sinA<p, which reduces 
to d 2 A(p w — uj 2 A(p from condition (c) that \A(p\ <C 1, 
where ui v — y/ J v /4irL\$>o\ 2 — v v / sqrtLS is the natural 
frequency of the given ISI junction. The solution in this 
case is 



Aip(t) ~ ipM s'm(u> v t + a) 



(14) 



with I^mI *C 1- In the presence of such vortex condensate 
phase oscillations, we may detect the magnetic vortex 
supercurrents by measuring the transverse voltage differ- 
ence across the region S (see Fig. |l|) which is induced by 
the vortex motion as follows: 
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2eAVtr = LJ X w LJ v ipM sin(w v i + a) 



(15) 



As the junction thickness <5 is increased above £ v , the 
Laplace equation of tp does not hold anymore in the re- 
gion S, and the relation J v cx 1/6 breaks down. There- 
fore, by varying the thickness 8 and studying J v (d), the 
correlation length £ v can be determined as the length 
scale of 8 where J v (8) deviates from J v oc 1/5. As we 
further underdope the charge-carriers of the region S, we 
expect that £ v will grow, resulting in a super-long-ranged 
dual proximity effect. 

We are now able to make estimates of various quanti- 
ties related to the (dual) proximity effects in the SIS (ISI) 
junctions, to see if we can find a regime of physical pa- 
rameters within experimental access. The information on 
one particular set of parameters is provided by the experi- 
mental results of Decca et al. || on the SIS junction. The 
insulating junction (with thickness S c ) described by Eq. 
@ leads to a Josephson current Ij = I c sin Atp where 
I c = en s L/2mi8 c , if 5 C < £ c . In Ref. §, I c w 2.6 fiA 
with S c ~ 45 nm and L w 200 nm, which leads to 
p s = n s /4mf w 1.2 meV. From this, using the ID Joseph- 
son weak link model ||i~9[| , one can estimate the magnetic 
vortex rest energy m v t> 2 = p s 7rmK « 2.6 x 10 -2 eV. 
Note that the above small magnitude of the superfluid 
phase stiffness p s already indicates that the material is 
under strong vortex fluctuations and the vortex mass m v 
may have been significantly renormalized since p s < 2 
meV ||, although we will keep the bare mass for order 
of magnitude estimates. 

Now we consider the ISI junction made of the same 
superconducting (S) and insulating (T and I2) regions 
as were used in Ref. Then the correlation length 
of the vortex fields in the region S is expected to be 
£ v sa (m v v v ) w 5.1 nm, which is a few times larger than 
the superconducting coherence length. We can enhance 
the correlation length £ v by further undoping the charge- 
carrier density. Here we assume a junction of length 
6 ~ 5.1 nm and explore the possibility of the simple so- 
lution in Eq. (|l4|). The magnitude of the fluctuations 
of $ in time [condition (a)] cannot be checked within 
the phenomenological theory presented here. Condition 
(b) (S < £ v and\dt®s/v F \ < l^^sl ) can be checked 
by comparing ti v |9a;$|/|$| ~ v v /S w 10 THz with the 
natural frequency |9 t $|/|$| ~ oj v = ^Jv^/L8 which is 
about 1 THz; we confirm that condition (b) can be sat- 
isfied with the given parameters. Condition (c) (|A<^| ~ 
AB/ (2(j) Lj v \&o\ 2 ) "C 1) gives the upper limit of the mag- 
netic field difference that allows the solution in Eq. (14). 
From £ c rs 100 nm and ^ 1 = 47r|$ | 2 w v ln(L/£ ), we 
find that the upper limit of AB is 2.2 x 10 2 Gauss. The 
upper limit of the transverse voltage in this case is about 
|AVtr| ~ LJ v ip M /2e - 11 mV, obtained from Eq. (|l|) 
and the relation (J v ^ c ) _1 = dln(L/£o)/v v . Therefore, 
based on the information obtained by Decca et al. 
the proposed experiment of dual Josephson junction and 



the observation of the solution in Eq. (|14j) may be within 
the reach of experimental access. 

We have shown that the anomalous proximity effect in 
an underdoped insulating cuprate Q can be interpreted 
in terms of vortex proliferation. We proposed that the 
dual counterpart of the proximity effect can be observed 
in an ISI junction experiment. We anticipate that a sim- 
ilar super-long-range proximity effect will be observed in 
the dual Josephson junction experiment as long as the 
junctions are prepared in close proximity to the SIT. 
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